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1 Introduction 


Evidence for the decay B° —> was first seen by CLEO H with a significance of 3.9 standard 

deviations (a). The decay was observed about a year ago with a sample of 89 million BB pairs 
by BABAR |2j with a significance of more than 7<r; the branching fraction was measured to be 
B(B° —► ujK°) = (5.9j]j;'3 ± 0.5) x 10 -6 . More recently Belle has published evidence (significance 
3.2a) for this decay mode with a sample of 85 million BB pairs [Hi, finding a branching fraction of 
(4.0^'g ± 0.5) x 10- 6 . Belle also has a preliminary study with a sample of 275 million BB pairs 
of the time-dependence of the decay with a sample of about 30 events [3j, though no branching 
fraction is reported. 

The world average branching fraction, (5.5lJ'f) x 10^ 6 [Sj, is somewhat larger than older theo¬ 
retical predictions El 13 El and a more recent prediction using QCD factorization E- A very recent 
paper m finds an enhancement of the QCD factorization prediction by more than a factor of two 
due to final-state interactions, in excellent agreement with the world average. A phenomenological 
fit that uses SU(3) flavor symmetry and all available measurements of pseudoscalar-vector (PV) 
decays (branching fraction and CP asymmetry measurements for more than 30 charmless decay 
modes) finds a branching fraction of ( 5 . 3 + 04 ) x 10 -6 for this decay jllj . 

In this paper we report improved branching fraction results for this decay as well as a measure¬ 
ment of the CP asymmetry parameters. In the Standard Model (SM), this decay is expected to 
proceed primarily through a penguin (loop) diagram as shown in Fig. |T|i. though a Cabibbo- and 
color-suppressed tree diagram is also possible (Fig. [TJj). Neglecting the suppressed amplitude, these 
decay modes have the same weak phase as the charmonium K° decays m which proceed through 
the Cabibbo-Kobayashi-Maskawa (CKM) favored b —> ccs amplitude. Thus the time-dependent 
asymmetry measurement for the decay B° —> ujK° would yield the same value of sin2 (3 as for the 
charmonium K° decays ITT] . Tests of this equality have been made from similar B° decays which 
are expected to be dominated by penguin amplitudes such as those to the charmless final states 
</>K°, rfK° , K + K~K°, ir°K° and / 0 (980)/i° JEJ. 



(a) (b) 

Figure 1: Feynman diagrams for the decay B° —► ujK°: (a) penguin diagram and (b) Cabibbo- and 
color-suppressed tree diagram. 


Additional higher-order amplitudes and non-SM amplitudes carrying different weak phases 
would lead to differences between the measurements of the time-dependent CP violating parame¬ 
ter in these rare decay modes and in the charmonium K° decays. The recent calculation involving 
final-state interactions that was mentioned above m predicts an increase in sin2/3 of about 0.10 for 
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B° —> uK° due to the color-suppressed amplitudes, though this increase is nullified when final-state 
interactions are included. 

2 The BABAR Detector and Dataset 

The results presented here are based on data collected with the BABAR detector [161 at the PEP-II 
asymmetric e + e~ collider located at the Stanford Linear Accelerator Center. We use a data sample 
with an integrated luminosity of 211 fb _1 recorded at the T(4S) resonance (center-of-mass energy 
y/s = 10.58 GeV). This corresponds to 232 million BB pairs. The asymmetric beam configuration 
in the laboratory frame provides a boost of j3' y = 0.56 to the T(45). 

Charged particles from the e + e - interactions are detected and their momenta measured by a 
combination of a vertex tracker (SVT) consisting of five layers of double-sided silicon microstrip 
detectors and a 40-layer central drift chamber, both operating in the 1.5-T magnetic field of a super¬ 
conducting solenoid. We identify photons and electrons using a CsI(Tl) electromagnetic calorime¬ 
ter (EMC). Further charged particle identification (PID) is provided by the average energy loss 
( dE/dx ) in the tracking devices and by an internally reflecting ring imaging Cherenkov detector 
(DIRC) covering the central region. The flux return of the solenoid is composed of multiple layers 
of iron and resistive plate chambers for the identification of muons and long-lived neutral hadrons. 

3 Time-dependent Analysis 

From a BB pair we reconstruct a B° or B° decaying into the CP eigenstate B° — > luK® (Bqp). We 
also reconstruct the vertex of the other B meson (Stag) and identify its flavor. The time difference 
At = tcp — ttag; where tcp and f tag are the proper decay times of the signal and tagged B mesons, 
respectively, is obtained from the measured distance between the Bcp and I3 tag decay vertices and 
from the boost (/Ty = 0.56) of the T(4S') system. The distribution of At without detector resolution 
effects is: 

e -\M\/r 

F(At) =---{1 Aw ± (1 — 2 w) [Ssm(ArridAt) — C cos(AmrfAt)]} , (1) 

4r 

where the upper (lower) sign denotes a decay accompanied by a B° (B°) tag, r is the B° lifetime 
Amd is the mixing frequency, and the rnistag parameters w and Aw are respectively the 
average and difference of the probabilities that a true B° (B°) meson is tagged as B° (B°). The 
tagging algorithm m has seven mutually exclusive tagging categories of differing purities (including 
one for untagged events that we retain for the branching fraction determination). Separate neural 
networks are trained to identify primary leptons, kaons, soft pions from D* decays, and high- 
momentum charged particles from B decays. Each event is assigned to one of these categories 
based on the estimated rnistag probability and on the source of tagging information. The measured 
analyzing power, equal to the reconstruction efficiency times (1 — 2 w) 2 summed over all categories, 
is (30.5 ± 0.6)%; this is determined from a large sample of L>-decays to fully reconstructed flavor 
eigenstates (-Bfl av ). The parameter C measures direct CP violation. If C = 0, then S = sin2/3 aside 
from the corrections discussed in the Introduction. 
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4 Event Selection and Analysis Method 


To establish the event selection criteria, we use Monte Carlo (MC) simulations [TS] of the signal 
decay modes, BB backgrounds, and detector response. We reconstruct B candidates by combining 
Kg and u ; candidates. We select K® —► decays by requiring the 7r + 7r~ invariant mass to be 

within 12 MeV of the nominal K® mass. We further require the three-dimensional flight distance 
from the beam spot to be greater than three times its uncertainty in a fit that requires consistency 
(fit probability greater than 0.001) between the flight and momentum directions. We reconstruct 
oj mesons through the primary oj —> 7r + 7r _ 7r° decay channel from two charged tracks and a 7r° 
candidate formed from pairs of photons with energy greater than 50 MeV and invariant mass 
between 120 and 150 MeV. The 7r + 7r _ 7r° invariant mass is required to be between 735 and 825 
MeV. For the time-dependent analysis, we require |At| < 20 ps and a At < 2.5 ps. We find an 
average of 1.13 B candidates per event. We choose the candidate with the 7r + 7r~7r° mass nearest 
to the nominal uj mass m- 

For a correctly reconstructed H-meson candidate, the mass must equal the nominal B mass and 
the reconstructed energy must be equal to one-half the center of mass energy. Thus we characterize 
a candidate kinematically by the energy-substituted mass mES = [(^s + Po ■ p b) 2 /Bo — p 2 B \i and 
energy difference A E = E* B — \\Js, where the subscripts 0 and B refer to the initial T(45) and 
to the B candidate, respectively, and the asterisk denotes the T(4S) rest frame. The resolution on 
A E (mEs) is about 30 MeV (3.0 MeV). We require \AE\ < 0.2 GeV and 5.25 < mES < 5.29 GeV, 
and include both of these observables in the maximum-likelihood (ML) fit (see Sec. El). 

5 Backgrounds 

To reject background from continuum e + e” —> qq events (q = u, d, s, c), we use the angle 6t between 
the thrust axis of the B candidate and that of the rest of the tracks and neutral clusters in the 
event, calculated in the center-of-mass frame. The distribution of cos 6t is sharply peaked near 
±1 for combinations drawn from jet-like qq pairs and is nearly uniform for the isotropic B meson 
decays; we require | cos 9t\ < 0.9. The remaining continuum background dominates the samples. 

We use MC simulations of B°B° and B + B~ production and decay to investigate BB back¬ 
grounds. We estimate that this background comprises 0.2% of the fit sample. Since we estimate 
from simulation studies that this background would change the signal yield by less than one event, 
we do not include a BB component in the fit. 

6 Maximum Likelihood Fit 

We use two unbinned, multivariate maximum-likelihood fits, one to extract signal yields and and 
one to determine the CP violating parameters. The yield fit does not use the tagging and At 
information in order to reduce systematic errors from the At parameterization (though the yields 
are in excellent agreement for both fits). We use six discriminating variables: m es, A E, At, the 
7T + 7r~7r° invariant mass (m^), a Fisher discriminant E, and TL = | cos0#j. The Fisher discriminant 
combines five variables: the polar angles, with respect to the beam axis in the T(45 l ) frame, of the 
B candidate momentum and of the B thrust axis; the tagging category; and the zeroth and second 
angular moments Lop of the energy flow about the B thrust axis. The moments are defined by 
Lfc = EiPi x l cos ®i\ k i where pi is the momentum of track or neutral cluster i, 0i is its angle in the 
T(45) frame with respect to the B thrust axis and the sum excludes the B candidate daughters. 
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The helicity angle 9h is the angle, in the u rest frame, between the normal to the u> decay plane 
and the B direction. For each species j (signal or background) and each tagging category c, we 
define a total probability density function (PDF) for event i as 

Vic = V j (mES i )Vj(AE i )V j (J ri )V j (rn i u) )V j ('H i )Vj(At i , a\ t , c) , (2) 

where a\ t is the error on At for the event i. With rij defined to be the number of events of species 
j and fj jC the fraction of events of species j for each category c, we write the extended likelihood 
function for all events belonging to category c as 

N c 

C c = exp ^ n jfj,cJ J^[(^sig/sig,c^ s ig,c + ^bkg/bkg,c7'’bkg )> (^) 

3 i 

where N c is the total number of input events in category c. The total likelihood function for all 
categories is given as the product over the tagging categories. For the yield-only fits, we integrate 
over the tagging categories and the product is over the total number of events in the sample. 

We maximize the likelihood function by varying a set of free parameters: S'; (7; signal and back¬ 
ground yields; background shape of F , AE, m^s, and vr + 7r _ 7r° mass; the fractions of background 
events in each tagging category; and six parameters representing the background At shape. We 
determine the PDF parameters for signal from simulation except for At, where we use the -Bfiav 
data sample discussed in Sec. E3 For the continuum background we use (mES; AE) sideband data 
to obtain initial values, before applying the fit to data in the signal region. We parameterize each 
of the functions T’sig(niEs)) 'Psig(A-Efc), Vj(F), ? s igK) and real ui component of PbkgK) with 
either a Gaussian, the sum of two Gaussians or an asymmetric Gaussian function as required to de¬ 
scribe the distribution. Slowly varying distributions (mass, energy or helicity-angle for continuum 
background and the combinatorial background component of Pbkg {mj)) are represented by linear 
or quadratic dependencies. The peaking and combinatorial components of the background 7r + 7r _ 7r° 
mass spectrum each have their own 7 ~i shapes. The continuum background in tues is described by 
the function x\Jl — x 2 exp [—£(1 — x 2 )\, with x = 2rriEs/\/s an d £ as a free parameter. The back¬ 
ground At shape is the sum of a core Gaussian convolved with an exponential function and two 
“tail” Gaussian functions. To verify the simulated resolutions in AE and tbeS; we use large control 
samples of the decays B~ —> tt~D° with D° —» which have a topology similar to that of 

the signal. Where the control data samples reveal differences from MC in tties ° r AE, we shift or 
scale the resolution function used in the likelihood fits. 

Before applying the fitting procedure to the data to extract the signal yields we subject it to 
several tests. Internal consistency is checked with fits to ensembles of “experiments” generated 
by MC from the PDFs. From these we establish the number of parameters associated with the 
PDF shapes that can be left free in addition to the yields. Ensemble distributions of the fitted 
parameters verify that the generated values are reproduced with the expected resolution. The 
ensemble distribution of ln£ itself provides a reference to check the goodness of fit of the final 
measurement once it has been performed. 

We evaluate biases from our neglect of correlations among discriminating variables in the PDFs 
by fitting ensembles of simulated experiments. Each simulated experiment has the same number 
of events as the data for both background and signal; background events are generated from the 
background PDFs while signal events are taken from the fully simulated MC samples. We find a 
positive bias of 7.3 ±0.5 events. Since events from a weighted mixture of simulated BB background 
decays are included, the bias we measure includes the effect of the neglect of BB background in 
the fit. 
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7 Results 


The results of the fits are shown in Table |T| The branching fraction is determined from the fit 
yield, corrected for the bias discussed above. The statistical error on the signal yield is equal to 
the change in value that corresponds to an increase of —21n£ by one unit from its minimum. 
The significance is equal to the square root of the difference between the value of —21n£ (with 
systematic uncertainties included) for zero signal and the value at its minimum. Results from 
simulated experiments suggest a possible underestimate of the fit error on S and C. To account 
for this effect, the statistical uncertainty for S and C have been increased by a factor of 1.07. 

Table 1: Results from yield and At fits. 


Quantity 

Yield fit 

At fit 

Events in ML fit 

9145 

8070 

Fit signal yield 

96 ± 14 

92 ± 13 

Efficiency (e) 

0.21 

- 

UBi 

0.31 

- 

e x n Bi 

0.065 

- 

Significance 

8.6 

- 

£>(10 -6 ) 

5.9 ± 1.0 ±0.4 

— 

S 

— 

o ^n +0 ' 34 

U - OU ^0.38 

C 

— 

—0 56 +0 ' 29 
U.OO_o.27 


In Fig. [3 we show projections onto tues and A E of a subset of the data for which the signal like¬ 
lihood (computed without the plotted variable) exceeds a threshold that optimizes the sensitivity. 
In Fig. El we show plots of the At distribution for B°- and R°-tagged events and their difference. 

8 Systematic Uncertainties and Crosschecks 

8.1 Branching Fraction Fit 

Most of the systematic uncertainties arising from lack of knowledge of the PDFs have been included 
in the statistical errors since most background parameters are free in the fits. For the signal the 
uncertainties in PDF parameters are estimated from the consistency of fits to MC and data in con¬ 
trol modes. Varying the signal PDF parameters within these errors, we estimate the uncertainties 
in the signal PDFs to be 0.7 events. The uncertainty in the fit bias correction is conservatively 
taken to be half of the correction itself. 

The above uncertainties are additive in nature. There are also multiplicative systematic er¬ 
rors that are comparable in size, primarily uncertainties in the efficiency. The latter, found from 
auxiliary studies, include an uncertainty in the absolute efficiencies for tracking (1.4%), 7r° recon¬ 
struction (3.0%), and K® reconstruction (2.1%). Our estimate of the systematic error in B counting 
is 1.1%. Published data PH provide the uncertainties in the w product branching fractions (1%). 
The uncertainty in the efficiency of the cos$t requirement is 0.5%. The total systematic error on 
the branching fraction is 0.4 x 10~ 6 . 
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Figure 2: (a) toes and (b) A E projections for B° —> uK® for data subsets optimized from the 
signal likelihood. Points with error bars represent the data, solid curves the full fit functions, and 
dashed curves the background functions. 


8.2 At Fit 

The contributions to the systematic uncertainties in S and C are summarized in Tabled We eval¬ 
uate the uncertainties associated with the PDF shapes by variation of the parameters describing 
each discriminating variable that is not free in the fit. Systematic errors associated with signal 
parameters (At resolution function, tagging fractions, and dilutions) are determined by varying 
their values within errors. Uncertainties due to Aand tb are obtained by varying these param¬ 
eters by the uncertainty in their world average values EH All changes are combined in quadrature 
obtaining an error of 0.01 for both S and C. 

We vary the SVT alignment parameters in the signal MC events by the size of misalignments 
found in the real data. The resulting uncertainty in both S and C is negligible. 

For some tag-side B decays, there is interference between the CKM-suppressed b —> ucd ampli¬ 
tude and the favored b —► cud amplitude. We use simulation, allowing the full range of variation of 
the relevant parameters, to estimate the systematic errors due to this effect to be negligible for S 
and 0.015 for C. An uncertainty of 0.02 in S and C is assigned to account for limitations of Monte 
Carlo statistics and modeling of the signal. The uncertainty in the effect of the neglect of the small 
BB background is estimated to be less than 0.01 for both S and C. We find that the effects of the 
uncertainty in the position and size of the beam spot are negligible. The total systematic error is 
obtained by summing individual errors in quadrature. 

When we fit with the value for C fixed to zero, we find a shift in S of 0.09, consistent with the 
correlation of ~20% between these variables. We produce samples of pseudo-experiments generated 
with events produced to match the PDF distributions. From these samples, we verify that the fit 
bias on S and C is negligible and that there is a good agreement between expected and observed 
errors. 
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Figure 3: Projections onto At, showing the data (points with errors), fit function (solid line), and 
background function (dashed line), for (a) B° and (b) B° tagged events and (c) the asymmetry 
between B° and B° tags. 
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Table 2: Estimates of systematic errors. 


Source of error 

a(S) 

a{C) 

PDF Shapes 

0.01 

0.01 

Tag-side interference 

0.00 

0.02 

At modeling 

0.02 

0.02 

BB background 

0.01 

0.01 

Total 

0.02 

0.03 


9 Conclusion 

We use a reconstructed signal sample of 96 B° —> ujK® events to determine the branching fraction, 
B(B° —> loK°) = (5.9 ± 1.0 ± 0.4) x 10~ 6 . This value is in good agreement with previous measure¬ 
ments and with the world-average value for the charged mode B(B + —► uK + ) = 5.1 ± 0.7 0 . This 
result is also in excellent agreement with the QCD factorization prediction modified by including 
final-state interactions m and with results from flavor-SU(3) fits to data for charmless B decays 
to PV final states EH- 

We also measure the time-dependent CP-violating parameters to be S = O.OO^q'^ ± 0.02 and 
C = — 0.56^0 27 ^ 0.03. The measurement of C is consistent with zero and S is in good agreement 
with the value of sin 2(5 as measured in charmonium K° decays |12j . The value for S is also in 
agreement with, but more precise than, the Belle measurement S = 0.75 ± 0.641q'^ g 0|. 
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